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T)  Abstract 


01:--  .V-c  ;.-r 


Considerable  progress  was  made  in  the  areas  of  application Of 


the  random-choice  method  to  explosion  and  implosion  dynamics  (flows 
with  changes  in  area  and  detonation)  and  weak  spherical  shock-wave 
transitions  (flows  with  viscosity,  heat  conduction,  and  vibrational 
excitation  -  a  first).  Some  progress  was  also  made  in  applying  the 
fluid-in-cell  technique  to  oblique-shock-wave  reflections  (OSWR). 

The  method  looks  promising.  A  great  deal  of  experimental  and 
analytical  work  was  achieved  on  OSWR  in  air  and  SFg.  The  dusty-gas 
shock  tube  is  being  assembled  at  a  slow  pace  owing  to  a  lack  of  funds 
(support  from  ARO  was  terminated  in  February  and  some  funds  promised 
by  AFWL  have  also  not  been  forthcoming.  Consequently,  Canadian  money 
that  could  have  gone  into  dusty-gas  shock-tube  equipment  is  now  being 
used  for  salaries  and  operating  expenses  of  a  reduced  research  and 
technical  staff).  Two  analytical  papers  on  dusty-gas  shock-tube  flows 
have  appeared  in  the  Proceedings  of  the  Royal  Society  of  London  with 
AFOSR  support,  as  well  as  other  papers  and  UTIAS  reports. 

2)  Objectives 

Our  objectives  were  to  provide  analytical,  numerical  and  experimental  da 
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OSWR  in  pure  and  dusty  gases.  The  results  can  prove  to  be  invaluable 
in  verifying  complex  two  and  three-dimensional  codes  used  for 
nuclear  explosions  and  height-of-burst  chemical-explosive  experiments. 
Our  isopycnic,  pressure  and  interferometric  data  can  establish  an 
important  base  for  testing  these  complex  codes  so  that  they  can 
become  credible  before  many  millions  of  dollars  are  spent  in  running 

routine  and  unproven  numerical  results. 

/ 

3.  Description  of  Accomplishments 

and 

4.  Publications 

No  better  elaboration  can  be  given  of  the  Abstract  on  the 
accomplishments  achieved  during  the  past  36  months  than  through  the 
various  publications.  A  summary  and  a  figure  of  each  is  enclosed: 

a)  "On  a  Dusty-Gas  Shock  Tube"  by  Miura  and  Glass,  Proc.  R.  Soc.  Lond., 
A382,  373-388,  198Z. 

b)  "On  the  Passage  of  a  Shock  Wave  Through  a  Dusty-Gas  Layer"  by 
Miura  and  Glass,  Proc.  R.  Soc.  Lond.,  A385,  85-105,  1983. 

c)  "Temperature  Measurements  at  an  Implosion  Focus"  by  Saito  and 
Glass,  Proc.  R.  Soc.  Lond.,  A384,  217-231,  1982. 

d)  "fin  Update  on  No.n-Stationary  Oblique  Shock-Wave  Reflections: 

Actual  Isopycnics  and  Numerical  Experiments"  by  Deschambault  and 
Glass,  J.  Fluid  Mech.,  to  be  published  1983. 

e)  "An  Assessment  of  Recent  Results  on  Pseudo-Stationary  Oblique- 

Shock-Wave  Reflections",  by  Shirouzu  and  Glass,  UTIAS  Report  No. 
264,  1982. 

f)  "Domains  and  Boundaries  of  Pseudo-Stationary  Oblique  Shock-Wave 
Reflections",  by  J.-H.  Lee  and  I.  I.  Glass,  UTIAS  Report  No. 

202,  1982. 

g)  "Random-Choice  Solutions  for  Spherical  Shock-Wave  Transitions 
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of  N-Waves  in  Air  with  Vibrational  Excitation"  by  H.  Honma  and 
I.  I.  Glass,  UTIAS  Report  No.  2b3,  1982. 

h)  "An  Experimental,  Analytical  and  Numerical  Study  of  Temperatures 
Near  Hemispherical  Implosion  Foci",  by  T.  Saito,  UTIAS  Report 
No.  260,  1982  -  Also  U.  of  T.  Thesis  (Ph.D.). 

1)  "Numer'r^?  Analysis  of  Dusty  Supersonic  Flow  Past  Remote 

Axisymmetric  Bodies",  by  H.  Sugayama,  UTIAS  Report  No.  267,  1983. 

5.  1  Professional  Personnel 

Prof.  I.  I.  Glass  -  Principal  Investigator 
Dr.  J.  J.  Gottlieb  -  Coinvestigator 
Prof.  H.  Honma  -  Visiting  Professor  (Japan) 

Dr.  J.-H.  Lee  -  Postdoctoral  Fellow  (Japan) 

Dr.  H.  Miura  -  Postdoctoral  Fellow  (Japan) 

Mr.  M.  Shirouzu  -  Senior  Research  Engineer  (Japan) 

Dr.  T.  Saito  -  Ph.D.  Candidate  and  Postdoctoral  Fellow  (Japan) 

M.A.Sc.  Candidates 
John  Hu  (Hong  Kong) 

John  Wheeler  (Canada) 

Fred  Wong  (Hong  Kong) 

Technicians  (Part  Time) 

Mr.  A.  Morte 
Mr.  D.  Wilmut 

6.  POD  Scientific  Interactions 

a)  Or.  Wes  Kitchens  and  Mr.  John  Keefer,  BRL,  can  verify  the  signific¬ 
ance  of  our  work  on  shock-wave  reflection  in  pure  and  dusty  gases  - 
especially  air,  in  providing  an  analytical,  numerical  and  experimental 
data  base. 

b)  The  same  can  be  said  of  Dr.  George  Ullrich,  DNA.Dr.  Wilson, 
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AFOSR,  and  Dr.  Harland  Glaz,  NSWC. 
c;  In  the  industrial  sector: Dr.  Allen  Kuhl  of  R  &  D  Associates, 

Dr.  Harold  Erode  of  Pacific  Sierra  Research,  Dr.  Charles  Needham, 
S-Cubed  and  Dr.  John  Kurylo,  Physics  International, 
d)  Additional  confirmation  can  be  obtained  from  Dr.  Phillip  Colella, 
Lawrence  Berkeley  Laboratory  and  Dr.  Paul  Kutler,  NASA  Ames. 

7.  Specific  Technical  Applications 

'  Without  our  work  unchecked  field-trial  explosion-data  and  numerical 
analyses  would  be  accepted  as  correct  at  face  value.  Now,  field- 
test  data  and  numerical  results  can  be  constructively  scrutinized 
by  comparison  to  our  work  and  the  correctness  of  the  physical  meaning 
of  all  data  can  be  properly  assessed. 


Proe.  R.  Sot.  Lend.  A  382.  373-333  (1932) 
Print'd  in  Great  Britain 


On  a  dusty-gas  shock  tube 

By  H.  MivRAf  and  1. I.  Glass 
Institute  for  Aerospace  Studies,  University  of  Toronto, 

4925  Dufferin  Street,  Doicnsview,  Ontario,  Canada  M3H  5T6 

( Communicated  by  A.  D.  Young,  F.R.S.  -  Received  10  December  1981) 

Analytical  and  numerical  methods  were  used  to  investigate  the  Row  in¬ 
duced  by  a  shock  wave  in  a  shock-tube  channel  containing  air  laden  with 
suspended  small  solid  particles.  Exact  results  are  given  for  the  frozen  and 
equilibrium  shock-wave  properties  as  a  function  of  diaphragm-pressure 
ratio  and  shock-wave  Mach  numbers.  The  driver  contained  air  at  high 
pressure.  A  modified  random-choice  method  together  with  an  operator- 
splitting  technique  show  clearly  both  the  decay  of  a  discontinuous  frozen 
shock  wave  and  a  contact  discontinuity,  and  the  formation  of  a  stationary 
shock  structure  and  an  effective  contact  front  of  finite  thickness. 

The  effects  of  particle  diameter,  particle-number  density  and  diaphragm- 
pressure  ratio  on  the  transitional  behaviour  of  the  flow  are  investigated  in 
detail.  The  alteration  of  the  flow  properties  owing  to  the  presence  of 
particles  is  discussed  thoroughly  and  compared  with  classical  shock -tube 
flows. 


1.  Introduction 

When  a  gas  carries  many  solid  particles,  they  significantly  affect  the  flow  through 
the  transfer  of  momentum  and  heat  from  or  to  the  gas.  Shock  waves  propagating 
in  such  a  dusty  gas  are  characterized  by  a  transition  region  orders  thicker  than 
that  caused  by  viscosity  and  heat  conduction  in  a  pure  gas.  Across  the  transition 
front,  the  interaction  of  the  gas  and  the  particles  leads  to  an  equilibrium  state  of 
the  mixture.  The  structure  of  stationary  shock  waves  has  been  studied  theoretically 
on  the  assumption  that  the  formulae  for  the  drag  and  the  rate  of  heat  transfer  for  a 
single  spherical  particle  placed  in  a  steady  flow  can  still  be  applied  to  the  motion 
of  many  particles  contained  in  a  dusty  gas  (Carrier  1958,  Marble  1963,  Kriebel 
1964,  Rudinger  1964,  Miura  1972). 

Shock  waves  were  studied  experimentally  in  a  dusty  gas  inside  a  shock  tube  to 
obtain  some  data  on  the  interaction  of  the  two  phases  (Crowe  el  al.  1963,  Selberg  & 
Nicholls  1968,  Rudinger  1970,  Outa  el  al.  1976).  Some  of  the  results  showed  an 
effective  drag  coefficient,  obtained  from  the  observation  of  the  acceleration  of  the 
particles  behind  the  shock  waves,  that  differed  appreciably  from  that  for  a  single 
particle.  However,  there  were  many  factors  influencing  the  results  and  a  definitive 
conclusion  on  the  appropriate  drag  coefficient  to  be  used  is  not  available  as  yet. 

Recently,  numerical  analyses  followed  these  experimental  studies.  The  shock- 

t  Permanent  address:  College  of  Engineering,  University  of  Osaka  Prefecture,  Sakai,  Japan. 
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On  a  dusty-gas  shock  tube 

lie  shock  front,  however,  are  still  fairly  large  at  this  stage.  These  must 
*  ith  time  (Miura  197a)  ami  the  extent  of  the  shock  transition  will  elongate. 
...t  the  shock  relaxation  length  or  time  is  larger  for  larger  mass  concentration. 

;•  mally,  the  effect  of  the  presence  of  particles  upon  the  stationary  shock  waves  is 
discussed  based  on  the  idealized  equilibrium -flow  limit.  The  pressure  pt  of  the  gas 
in  the  uniform  region  behind  the  shock-transition  region  satisfies  the  classical 
shock-tube  equation  (Glass  &  Hall  1959) 

gjg.L  l-****-”  (25) 

Pi  pA  (2y,)*{2y1  +  (y»+  - 


Figure  5.  Flow  quantities  at  T  *=  32  (a  —  I,  P„  =  10,  d  ■=  40  um):  (a)  temperature 
(6)  velocity.  Thick  lines,  gas;  thin  lines,  particles;  dashed  lines,  equilibrium  flow. 


Figure  6.  Flow  quantities  at  r  =  4  (*  *•  2,  P4l  «•  10,d  ■  10  pm):  (a)  temperature 
lb)  velocity.  Thick  lines,  gas;  thin  lines,  particles;  dashed  lines,  frozen  flow. 


Free.  B.  Soe.  Loud.  A  385,  85-105  (1983) 
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On  the  passage  of  a  shock  wave  through  a  dusty-gas  layer 

By  H.  MiuRAf  and  1. 1.  Glass}: 
f  College  of  Engineering,  University  of  Osaka  Prefecture,  Sakai,  Japan 
t  Institute  for  Aerospace  Studies,  University  of  Toronto,  Toronto,  Canada 

(Communicated  by  A.'D.  Young,  F.RJS.- Received  20  May  1982) 

[Plato  1] 

■  /  '  The  flow  resulting  from  the  passage  of  a  shock  wave  through  a  dusty-gas 
layer  is  studied  theoret  ically.  On  the  basis  of  an  idealized  equilibrium-gas 
approximation,  the  criteria  for  the  wave  reflexion  at  the  contact  surface 
separating  the  pure  gas  from  the  dusty-gas  layer  are  obtained  in  terms  of 
the  properties  of  the  gas  and  the  dusty  gas.  For  the  cases  treated  here,  a 
shock  wave  is  reflected  at  the  first  contact  surface  and  a  shock  wave 
stronger  than  the  incident  one  is  transmitted  into  the  dusty-air  layer. 
Subsequently,  a  rarefaction  wave  is  reflected  at  the  second  contact  surface 
and  the  shock  wave  transmitted  into  the  free  air  is  weakened  by  this  non¬ 
linear  interaction.  The  induced  rarefaction  wave  reflects  later  at  the  first 
contact  surface  as  a  compression  wave,  which  runs  through  the  layer  to 
overtake  the  transmitted  shock  wave  in  air.  The  final  emergent  shock 
wave  from  the  dusty  air  has  almost  the  same  strength  as  the  original 
shock  wave  entering  the  layer. 

The  time-dependent  transition  properties  through  the  shock  wave3, 
contact  surfaces  and  rarefaction  waves  are  found  by  solving  the  equations 
of  motion  numerically  by  a  modified  random-choice  method  with  an 
operator-splitting  technique. 

1.  Introduction 

High  speed  flows  of  a  mixture  of  a  gas  and  small  solid  particles  are  encountered  in 
several  brandies  of  engineering  and  sdence  (see,  for  example.  Marble  1970, 
Boothroyd  1971 ;  Rudinger  1980).  Shock  waves  in  such  a  dusty  gas  exhibit  relaxa¬ 
tion  features  typical  of  two-phase  flow.  For  sufficiently  strong  shock  waves,  a 
discontinuous  jump  in  the  gas  phase  precedes  a  thick  transition  region.  A  non¬ 
equilibrium  state  of  the  mixture  gives  way  gradually  to  an  equilibrium  transition 
through  the  transfer  of  momentum  and  heat  between  the  gas  and  the  partides. 
This  feature  was  used  in  experimental  studies,  with  a  shock  tube,  for  an  under¬ 
standing  of  the  basic  mechanisms  of  the  interaction  between  the  two  phases 
(Crowe  el  al.  1963,  Selberg  &  Nicholls  1968,  Rudinger  1970).  For  example,  effective 
drag  coefficients  of  the  particles  were  deduced  from  the  experimental  results,  but 
their  values  cannot  as  yet  be  considered  definitive. 

Some  experiments  used  a  shock  tube  in  which  only  a  limited  port  ion  of  the  channel 
away  from  a  diaphragm  was  filled  with  a  dusty  gas  (Outa  et  al.  1976,  Lowenstein  & 
♦  [  W]  VoL jsj.  A 
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Temperature  measurements  at  an  implosion  focus 
By  T.  Saito  asd  1. 1.  Glass 
Institute  for  Aerospace  Studies,  University  of  Toronto, 

■  4925  Dufferin  Street,  Doivnsviexr,  Ontario,  Canada  3I3H  5T6 

( Communicated  by  A.  G.  Guydon,  F.R.S.  -  Received  28  May  1982) 

//  * 

Spectroscopic  temperature  measurements  were  made  at  the  focal  point  of 
imploding  shock  waves  in  the  UTIAS  implosion  chamber,  which  has  a 
20cm  diameter  hemispherical  cavity.  The  chamber  was  filled  with  a 
stoichiometric  H2-08  gas  mixture  at  different  initial  pressures  {1.4— 

6.9  MPa).  The  mixture  was  ignited  at  the  origin  by  an  exploding  wire 
generating  an  outgoing  detonation  wave,  which  reflected  at  the  chamber 
wall  as  an  imploding  shock  wave  (gas  runs).  Experiments  with  an  explosive 
shell  of  pentaerythritol  tetranitrate  (PETN  :C5H8N4Oh)  placed  at  the 
hemispherical  wall  were  also  conducted.  The  shell  was  detonated  by  the 
impact  of  the  reflected  gaseous  detonation  wave  at  its  surface,  an  intense 
implosion  wave  being  generated  thereby  (explosive  run).  The  temperatures 
were  measured  at  the  implosion  focus  by  using  a  medium  quartz  Hilger 
spectrograph  with  an  eight-photocell  polvchromator  attachment  over  the 
visible  wavelength  range.  The  measured  radiation  intensity  distributions 
were  fitted  to  black -body  curves.  The  temperatures  were  1 6  000-13000  K, 
for  gas  runs,  and  15000-17000K,  for  explosive  runs.  The  continuous 
spectra  from  photographic  film  and  the  measured  emissivities,  which  were 
very  close  to  unity,  confirmed  that  the  plasma  was  a  black  body.  Numerical 
studies  with  the  random  choice  method  and  classical  strong-shock  theory 
were  used  to  analyse  the  flows  in  the  entire  range  of  the  implosion  process. 
They  provided  much  information  on  theentire  implosion  process  within  the 
restriction  of  a  perfect-gas  assumption,  which  was  found  to  be  reasonable 
as  a  first  step  in  this  kind  of  analysis.  The  experimental  data  were  com¬ 
pared  with  the  analytical  results.  For  both  gas  and  explosive  runs,  the 
temperatures  were  lower  than  the  calculated  values  and  reasonable 
explanations  are  given  for  this  deviation. 

Introduction 

The  UTIAS  explosive-driven  implosion  chamber  has  been  used  in  several  research 
areas  since  it  was  conceived  by  Class  (1972)  in  the  early  1900$,  The  extremely  high 
pressures  and  temperatures  generated  at  the  focus  of  implosions  were  used  for 
driving  projectiles  to  hypervelocities  (Flagg  k  Class  1968),  as  a  shock-tube  driver 
(Glass  et  al.  1974),  and  to  produce  diamonds  (Glass  k  Sharma  1976).  Most  recently, 
experimental  studies  were  done  on  deuterium-deuterium  fusion  reactions  ancl  on 
measurements  of  the  neutrons  and  y-rays  generated  by  the  reaction  (Glass  k  Sagie 
1982). 

The  conditions  produced  by  imploding  spherical  shock  waves  were  first  obtained 
by  Guderley  (1942).  His  self-similar  solutions,  however,  deal  only  with  very  strong 
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in  y  is  very  difficult  to  take  into  account  and  in  the  present  study  it  was  assumed  that 
y  was  constant  (y  =  1.14)  in  the  region  of  interest  (ca.  0.55  mm).  Real-gas  effects  at 
this  point  were  shown  to  have  a  negligible  effect. 

Comparison  of  experimental  results  and  analysis 
Figure  12  compares  the  experimental  results  and  the  expected  average  tem¬ 
peratures  over  a  0.55  mm  diameter  circular  area  at  the  origin.  The  calculated 
temperature  depends  on  the  position  of  the  imploding  shock  wave.  When  the  wave 
*'  radius  is  larger  than  the  radius  of  the  observation  area  only  the  wavefront  can  be 
seen  through  the  window.  When  the  imploding  shock  wave  comes  inside  the  obser¬ 
vation  area,  not  only  the  wavefront  but  also  the  region  behind  the  imploding  shock 
wave  is  seen  and  the  average  temperature  becomes  lower  than  the  temperature  at 


p  ,/MPa  PETN  shell  thickness/nun 


Figure  12.  Comparison  of  experimental  results  with  analysis  for 
(o)  gas  runs  and  (6)  explosive  runs. 

the  wavefront.  This  arises  from  the  rapidly  falli’ng  temperature  behind  the  implosion 
wave.  As  a  result,  the  average  temperature  has  a  maximum  when  the  radius  of  the 
imploding  shock  wave  is  slightly  smaller  than  the  radius  of  the  observation  area. 
The  radius  of  the  imploding  shock  at  which  this  occurs  is  called  l\*,  as  shown  in 
figure  12.  In  other  words,  there  is  a  certain  maximum  observable  temperature  for 
a  certain  size  of  observation  area  no  matter  how  small  the  imploding  shock  ware 
becomes.  This  maximum  temperature  and  the  temperatures  expected  to  be  observed 
when  the  imploding  waves  are  at  different  radii  from  the  centre  are  also  shown  in 
figure  12.  In  principle,  if  it  was  possible  to  place  the  observation  area  exactly  at  the 
focus  of  the  implosion  then  higher  temperatures  would  be  measured  with  decreasing 
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Nonstationary  oblique  shock-wave  reflections  over  compressive  wedges  in  air  and 
argon  were  investigated  using  infinite-fringe  interferometric  techniques.  These 
allowed  direct,  continuous  and  accurate  observations  of  the  isopycnics  (lines  of 
constant  density)  of  the  flow  field.  The  initial  pressures  for  these  experiments  were 
made  as  high  as  possible  (30  to  250  torr)  in  order  to  increase  the  number  of  isopycnics 
and  to  enhance  their  details  and  distribution  along  the  wedge  surface  over  a 
shock-Mach-number  range  2.0  <  J/s  $  8.7.  Included  in  the  study  were  two  cases  of 
regular  reflection  (RR)  and  one  of  each  of  single-Mach  reflection  (SMR),  complex-Mach 
reflection  (CMR)  and  double-Mach  reflection  (DMR)  for  air,  and  one  RR,  SMR,  OMR 
and  DMR  for  argon.  These  particular  cases,  which  we  investigated  previously  in  N, 
and  At  using  a  finite-fringe  technique,  have  been  used  by  computational  fluid 
dvnamicists  to  check  their  finite-difference  results  against  our  experimental  data.  It 
will  be  shown  that  the  isopycnic  structure  previously  reported  by  us  differs  in  detail, 
in  most  cases,  from  that  of  the  present  study.  The  major  difference  arises  from  the 
fact  that  it  was  only  possible  previously  to  obtain  discrete  points  on  isopycnics  and 
along  the  wedge  surface.  Consequently,  the  results  obtained  before  were  not  as 
accurate.  Comparisons  were  made  of  actual  wall-density  distributions  with  numerical 
simulations  of  the  density  contours  of  the  various  flows  obtained  by  a  number  of 
authors.  Each  numerical  method  displays  its  advantages  and  disadvantages  in 
describing  the  details  of  the  flow  fields.  The  present  experimental  results  for  air  are 
new.  They  are  of  great  interest  from  a  practical  viewpoint.  The  experiments  in  argon 
were  redone  to  provide  better  data  for  a  gas  free  from  real-gas  effects  in  the  range 
of  initial  conditions  considered ,  in  order  to  simplify  the  computations  in  the  numerical 
simulations.  Although  the  recent  numerical  simulations  are  better  than  those 
reported  previously,  additional  efforts  are  required  to  improve  the  predictions  of  the 
shape,  location  and  values  of  the  isopycnics  and  other  flow  isolines  in  the  various 
regions  and  along  the  wall,  and  to  render  the  predictions  free  of  computer  ‘noise*. 
It  is  worth  noting  that  real-gas  effects  did  not  play  any  role  in  determining  the  various 
transition  lines  of  RR,  SMR,  CMR  and  DMR;  a  different  claim  was  made  in  our 
previous  work.  Relaxation  of  nitrogen  in  air  can  be  seen  at  the  highest  shock  Mach 
numbers  (J/t  =  7.19  and  8.70),  with  relaxation  lengths  in  good  agreement  with 
accepted  predictions. 
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FroiTRE  2.  Isopycnic*  for  regular  reflection  (case  l):  (o)  present  experimental  result*,  if,  »»2. 05. 
dw  **  63.4°.  p%  *  230  tore,  p%  *  3.8"  x  10"‘  g/cm’,  T,  *■  298.4  K,  A  =*  6943  A;  (3)  Ando's  (1981) 
experimental  results  in  nitrogen,  if,  =  t.97,  =  63.4°,  p,  =>  50  tore,  r,  «»  SWi  K, 

pt  *»  7.58  x  10'*  g/cm*,A  =  3472  A;  (e)isopvcnics  obtained  from  fringe-shift  analysis  of  finite-fringe 
interferograma  by  Ben-Dor  &  Glass  (1978),  .V,  —  2.01 .  0W  «»  63. 4*.  p,  ««  30  tore.  7i  **  298.6  K;  (d\ 
numerical  simulation  of  the  contour*  of  the  density  Held  by  Schnever  (1975);  (e)  numerical 
simulation  of  the  contour*  of  the  density  field  by  Kutler  &  Shankar  (1977). 
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Summary 


the  assumptions  and  criteria  used  in  existing  analyses  in  determining  the 
regions  and  transition  lines  of  pseudo-stationary  oblique-shock -wave  reflections 
have  been  re-examined  in  order  to  improve  the  agreement  between  experiments  and 
computed  data  for  regular  (RR) ,  single-Mach  (SMR) ,  coropl'ex-Mach  (CMR)  and  double- 
Mach  reflection  (OMR}. 

It  is  shown  that  the  relaxation  lengths  for  vibration  and  dissociation  deter- 
mine  whether  frozen  or  equilibrium  gas  transition  lines  are  applicable.  For  example, 
at  an  initial  temperature  of  300  K  and  a  pressure  of  15  torr  (where  much  previams 
work  was  done)  ar.  '•  "h-M.,. iysis  would  not  be  required  for  shock  Mach 
numbers  Ms  <  9  ,  <  8  in  air  and  Ms  <  3  in  CO,. 

Yet,  the  -  ..:i  data  in  N, ,  CO?  and  very  recent  results  for 

air,  which  ar.  •  the  .lerion  (consistent  with  relaxation  lengths)  of  the 

angle  6,  bet*.,: . :  incident  and  reflected  shock  wave,  do  not  conclusively  support 

the  frozen  or  librium  gas  calculations  for  N,  and  air.  It  does  support  CO,  as 
an  equilibrium  gas  contrary  to  a  previous  conclusion  of  agreement  with  y  »  1.29. 

A  new  additional  and  necessary  criterion  for  the  transition  from  single  to 
complex  Mach  reflection  improves  the  agreement  between  analysis  and  experiment  amd 
is  consistent  with  the  requirements  of  the  relaxation  length  and  the  angle  6. 
However,  it  now  appears  that  a  more  accurate  criterion  is  required  for  the  boundary 
line  between  CMR  and  OMR.  •' 

A  more  detailed  examination  of  the  boundary-layer-displacement  slope  at  the 
point  of  regular  reflection  appears  to  eliminate  the  so-called  von  Neumann  paradox, 
and  explains  the  persistence  of  regular  reflection  below  the  transition  line  for 
the  occurrence  of  Mach  reflection. 

It  is  also  shown  that  at  the  triple  point  the  Mach  stem  can  vary  from  being 
perpendicular  to  the  wedge  surface  in  actual  experiments  by  as  much  as  -3.0*  to 
7.5*.  Consequently,  calculations  of  the  triple-point-trajectory  angle  x  °n  the 
basis  that  the  stem  is  perpendicular  is  not  always  well  founded. 

t 

It  is  verified  that  at  lower  shock  Mach  numbers  Ms  and  large  wedge  angles  %, 
the  experimental  evidence  shows  that  the  transition  lines  for  SMR  *  CMR  and 
CMR  *  DMR  converge  at  a  point  on  the  RR  *  MR  line,  contrary  to  a  previous  simplified 
analysis. 


SHOCK  MACH  NUMBER  Ms 

VARIATION  OF  ANGLE  6  VS  SHOCK-MACH  NUMBER  AT  FIXED  0^  FOR  Ar 

-  FROZEN  Ar  (y  «  1.667),  +  EXPERIMENTAL  DATA  FOR  6^ 

(REF.  4). 
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Summary 


The  reflection  of  oblique  shock  waves  in  air  in  pseudo-stationary 
flow  was  investigated  analytically  and  numerically.  The  transition 
boundaries  between  the  four  types  of  shock-wave  reflection  (regular 
RR,  single  Mach  SMR,  complex  Mach  CMR  and  double  Mach  DMR)  were 
established  up  to  M$  *  20  for  both  perfect  and  imperfect  air  in 
thermodynamic  equilibrium  (rotation-vibration  coupling,  vibrational 
excitation,  dissociation,  electronic  excitation  and  ionization). 

In  addition,  an  analysis  was  made  for  perfect  gases  with  differing 
Y,  in  order  to  clarify  the  effects  of  the  specific-heat  ratio  Y»  on 
the  shock -wave  configurations.  It  was  verified  that  the  reflected 
wave  angle  u'  was  a  very  sensitive  function  of  y,  and  a  decrease  in 
Y  lowered  the  value  of  u'  significantly  and  even  shifted  the  value 
of  w*  towards  negative  values  under  certain  conditions  of  Mach 
reflection.  This  phenomenon  occurred  in  a  perfect  gas  with  y  less 
than  1.4  and  in  imperfect  air.  However,  it  was  absent  in  perfect 
air  with  y  *  1.4.  *-•  'ton  of  the  perfect  and  imperfect  air 

results  shows  -  .  transition  lines  and  removes  the  con-  . 

jecture  of  •  t. ^.e-Mach  reflection. 

The  present  analytical  results  were  compared  with  the  available 
experimental  data  for  air  and  nitrogen  for  shock  Mach  numbers  up  to 
10.  From  the  available  experimental  cases,  it  was  clarified,  by 
examining  the  relaxation  lengths  behind  the  shock  waves,  that  the 
flow  states  behind  the  shock  fronts,  which  determined  the  wave 
systems  for  the  shock-wave  reflection,  were  frozen  or  nearly  frozen 
regarding  vibrational  excitation  and  dissociation.  Consequently, 
in  general,  the  present  perfect-gas  analysis  agreed  with  experiment. 
However  RR  persisted,  to  some  extent,  below  the  perfect-gas  termina¬ 
tion  line  determined  by  the  detachment  criterion,  and  SMR  and  DMR 
did  sometimes  occur  outside  their  analytically  predicted  domains. 

The  development  of  more  accurate  transition  criteria,  improvements  ’ 
in  the  methods  of  predicting  the  first  triple-point  trajectory 
angle,  accurate  locations  of  the  kink  and  triple  points  along  with 
possible. boundary-layer  effects  on  the  persistence  of  regular  1 
reflection  are  problems  to  be  resolved  in  the  future. 
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FIG.  23  COMPARISON  OF  PREDICTED  DOMAINS  OF  VARIOUS  TYPES  OF  REFLECTION 
AND  THE  EXPERIMENT  IN  THE  (Ms-6w)  PLANE.  SOLID  LINES  ARE  FOR 
PERFECT  AIR  WITH  y  =  1.40.  DASHED  LINES  ARE  FOR  IMPERFECT  AIR 
AT  P0  -  15  TORR  AND  T0  *  300  K.  SOLID  SYMBOLS  ARE  AIR  DATA 
FROM  DESCHAMBAULT  [REF.  39],  OPEN  SYMBOLS  ARE  N,  DATA  FROM 
BEN-DOR  [REF.  403.  RR  DU  ,  SMR  a  a  ,  CMR  0  0,  DMR  0  +. 
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FIG.  24  COMPARISON  OF  PREDICTED  DOMAINS  OF  VARIOUS  TYPES  OF  REFLECTION  AND 
THE  EXPERIMENT  TN  THE  (Ms-ew)  PLANE  IN  THE  RANGE  1.0  <  Ms  S  3.0.  . 

SOLID  LINES - ARE  FOR  PERFECT  AIR  WITH  y  =  1 .40.  DASHED  LINES 

-  ARE  FOR  IMPERFECT  AIR  AT  P0  a  15  TORR  AND  T0  -  300  K.  DASH-DOT 

LINE - SHOWS  PERSISTENCE  OF  REGULAR  REFLECTION  RR.  EXPERIMENTAL 

DATA  -  AIR:  SMITH  0  O  [REF.  9],  BAZHENOVA  ET  AL  D  [REF.  18], 
DESCHAMBAULT  OiTD  [REF.  39];  NITROGEN:  BAZHENOVA  ET  AL  A  [REF.  18] 
BEN-DOR  OAV  [REF.  40];  OXYGEN:  LAW  AND  GLASS  AY  [REF.  19].  RR  t 
SMR  ▲  A  A  A  ,  CMR  ▼  v  ▼  ,  DMR  n  D  ,  MR 


RANDOM-CHOICE  SOLUTIONS 

FOR  WEAK  SPHERICAL  SHOCK-WAVE  TRANSITIONS  OF  K- WAVES  IN  AIR 
WITH  VIBRATIONAL  EXCITATION 

by 

H.  Honma  and  I.  I.  Glass 

Summary 


In  order  to  clarify  the  effects  of  vibrational  excitation  on  shock-wave  transi¬ 
tions  of  weak,  spherical  N-waves,  which  were  generated  by  using  sparks  and  exploding 
wires  as  sources,  the  compressible  Navier-Stokes  equations  were  solved  numerically, 
including  a  one-mode  vibrational -relaxation  equation.  A  small  pressurized  air-sphere 
explosion  was  used  to  simulate  the  N-waves  generated  from  the  actual  sources.  By 
employing  the  random-choice  method  (RCM)  with  an  operator-splitting  technique,  the 
effects  of  artificial  viscosity  appearing  in  finite-difference  schemes  were  eliminated 
and  accurate  profiles  of  the  shock  transitions  were  obtained.  However,  a  slight  ran¬ 
domness  in  the  variation  of  the  shock  thickness  retrains.  It  is  shown  that  a  computer 
simulation  is  possible  by  using  a  proper  choice  of  initial  parameters  to  obtain  the 
variations  of  the  N-wave  overpressure  and  half-duration  with  distance  from  the  source. 
The  calculated  rise  times  are  also  shown  to  simulate  both  spark  and  exploding-wire 
data.  It  was  found  that,  in  addition  to  the  vibrational-relaxation  time  of  oxygen', 
both  the  duration  and  the  attenuation  rate  of  a  spherical  N-wave  are  important  factors 
controlling  its  rise  time. 

The  effects  of  the  duration  and  the  attenuation  rate  of  a  spherical  N-wave  on  its 
rise  time,  which  are  designated  as  the  ll-uave  effect  and  the  nonstatiemary  effect, 
respectively,  are  discussed  in  more  detail  pertaining  to  Lighthill’s  analytical  solu¬ 
tions  and  the  RCM  solutions  for  nonstationary  plane  waves  and  spherical  N-waves.  It 
is  also  shown  that  the  duration  and  the  attenuation  rate  of  a  spherical  N-wave  are 
affected  by  viscosity  and  vibrational  nonequilibrium,  so  that  they  can  deviate  from 
the  results  of  classical,  linear  acoustic  theory  for  very  weak  spherical  waves. 
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SUMMARY 


Spectroscopic  temperature  measurements  in  the  visible  radiation 
range  vere  made  near  the  focus  of  imploding  shock  waves  in  the  UTIAS 
implosion  chamber,  which  has  a  20-cm  diameter  hemispherical  cavity.  The 
chamber  was  filled  with  a  stoichiometric  Hg-Og  gas  mixture  at  different 
initial  pressures  (13.6  -  68.1  atm).  The  mixture  was  ignited  at  the  origin 
by  an  exploding  wi  -rating  an  outgoing  detonation  wave  which  reflected 

at  the  chm  un  imploding  shock  wave  (gas-runs).  Additional 

experice-!: .  -nere  an  explosive  shell  of  PETN  was  placed  at  the  hemispherical 
wall  vere  also  conducted.  The  shell  was  detonated  by  the  impact  of  the 
reflected  gaseous  detonation  wave  at  its  surface,  thereby  generating  an 

.  -*1 

intense  implosion  wave  (explosive-runs).  The  continuous  spectra  and  the 
measured  emissivities,  which  vere  close  to  unity,  confirmed  the  blackbody 

i 

nature  of  the  gas  near  the  implosion  focus.  Temperatures  were  determined 

•  . 
by  fitting  the  measured  radiation-intensity  distributions  to  Planck's 

radiation  distribution  curve.  They  were  10,000  -  13,000  K  for  toe  gas 

and  15,000  *  17,000  K  for  explosive  runs. 

Semi-analytical  solutions  vere  developed  by  combining  the  solutions 
of  the  random-choice  method  (ROM)  and  the  classical  strong-shock  theory 
of  Guderley.  The  semi-analytical  solutions  provided  a  great  deal  of  infor¬ 
mation  on  the  entire  implosion  .process  within  the  restriction  of  a  perfect- 
gas  assumption  which  was  found  to  be  reasonable  as  a  first  step  in  this 
kind  of  analysis.  The  RCM  alone  detailed  the  ignition  process  of  a  PETN 
explosive  shell  at  the  hemispherical  implosion-chamber  vail  and  the  sub¬ 
sequent  complex  wave  interactions  in  that  region. 

A  comparison  between  the  analysis  and  the  experiments  showed  that 
the  observed  temperatures  for  the  explosive  runs  vere  limited  to  relativelj 


iii 


low  values  (15,000  ~  17,000  K)  owing  to  the  screening  effect  of  the 
preheated  gas  layers  ahead  of  the  radiating  imploding  shock  waves, 
they  converged  to  very  small  radii. 


rf»- ’ 


NUMERICAL  ANALYSIS  OF  DUSTY  SUPERSONIC  FLOW 
PAST  BLUNT  AXISYMMETRIC  BODIES  „ 


by 


H.  SUGIYAMA 


Submitted  Deceaber,  1982, 

» 


March,  1983 


UTIAS  Report  No.  267 
CN  ISSN  0082-S25S 


Summary 


An  inverse  method  was  developed  for  treating  gas-particle  supersonic  flow 
past  axi symmetric  blunt  bodies.  This  method  is  based  on  two  transformations 
(von  Xises  and  an  additional  one),  which  are  convenient  for  determining  the  • 
shock-layer  flow  fields  and  the  body  shapes. 

In  using  the  present  method,  the  pure  gas  flow  fields  around  spheres  were 
first  solved  numerically  for  the  freestream  Mach  numbers  ■  10,  6,  4,  5,  2 
and  l.S.  Tnese  were  found  to  be  in  very  good  agreement  with  the  available 
results  of  Van  Dyke  and  Gordon.  Then  the  gas-solid-particle  flow  in  the  shock 
layer  around  blunt  bodies  (nearly  spheres)  were  solved  for  the  freestream  Mach 
numbers  ■  10  and  l.S,  with  freestream  loading  ratios  o  ■  0,  0.2,  0.5  and 
1,0  and  particle  diameters  dp  ■  1,  2,  5  end  10  vm,  respectively.  The  effects 
of  M^,  2„  end  o  on  the  shock-layer  thickness  and  the  body  surface  pressures  are 
discussed.  The  variations  of  the  flow  properties  along  the  stagnation  and 
adjacent  streamlines  are  also  shown  in  some  detail. 


